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SUMMARY 


A  three  dimensional  finite  element  analysis  of  cracks  in  the  MB326H  Macchi  wing  lower  spar 
c<^>  is  presented.  The  case  of  cracks  emanating  from  the  intersection  of  a  wing  attachment 
fitting  bolt  hole  and  a  web  attachment  fastener  hole  at  wing  station  Y900  is  considered. 
Stress  intensity  factors  for  two  proposed  crack  geometries  and  the  residual  strength  of  the 
ligament  above  the  web  attachment  fastener  hole  in  the  spar  cap  flange  are  calculated. 
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1. 


INTRODUCTION 


The  wing  lower  spar  cap  of  the  Macchi  MB326H  is  susceptible  to  fatigue  cracking 
at  wing  station  Y910,  near  the  end  of  the  wing  attachment  fitting  (WAF).  A  life 
cissessment  of  the  spar  cap  was  completed  in  April  1992  [1].  This  was  concerned 
primarily  with  cracks  emanating  from  a  web  attachment  fastener  hole  located  in  the 
flange  of  the  spar  cap. 

Following  the  repair  of  several  wing  lower  spar  caps  it  became  apparent  that  some 
of  the  blind  web  attachment  holes  in  the  base  of  the  spar  cap  ended  in  close  prox¬ 
imity  to  WAF  bolt  holes  and  that  this  may  influence  the  fatigue  life  of  the  spars. 
Although  no  crack  growth  from  these  nearly-intersecting  holes  has  been  observed, 
it  was  considered  that  further  investigation  into  the  problem  Wcis  appropriate  given 
the  history  of  fatigue  cracking  in  the  spar  cap  and  the  difficulty  of  inspection.  An 
alternative  crack  growth  scenario  was  postulated,  see  [2],  which  involved  cracks  ini¬ 
tiating  from  the  intersection  of  a  WAF  bolt  hole  and  a  web  attachment  fastener 
hole.  It  was  also  postulated  that  such  cracking  may  cause  failure  in  the  ligament 
above  the  web  attachment  fastener  hole  in  the  spar  cap  flange  and  that  this  may  be 
used  as  a  basis  for  inspection. 

It  was  recognised,  however,  that  it  may  also  be  possible  for  cracks  growing  from 
intersecting  holes  to  progress  to  failure  of  the  section  prior  to  failure  of  the  adj:'.cent 
uncracked  flange  hole  ligament.  In  this  case,  flange  surface  NDI  would  detect  no 
flange  hole  cracking  and  would  provide  no  assurance  of  structural  safety.  To  check 
this  possibility,  two  intersecting  hole  crack  cases  (with  no  spar  flange  cracking)  were 
investigated,  the  residual  strengths  of  the  cracked  regions  being  compared  with 
that  of  the  (intact)  flange  hole  ligament.  The  site  likely  to  fail  first  could  thus  be 
ascertained. 

A  three-dimensional  finite  element  analysis  of  the  lower  spar  cap  and  WAF  assembly 
between  wing  stations  Y884  and  Y920  was  undertaken.  The  residual  strength  was 
determined  for  the  uncracked  structure  and  for  the  structure  containing  assumed 
crack  geometries  of  two  sizes. 


2.  FINITE  ELEMENT  MODEL 

The  end  of  the  lower  spar  cap  and  wing  attachment  fitting  assembly  between  wing 
stations  Y884  and  Y920  was  modelled.  This  included  the  last  three  WAF  bolts, 
see  Figure  1.  The  web  attachment  fastener  holes  (Deutsch  fasteners)  D15  and  B16 
are  adjacent  to  WAF  bolt  hole  C7.  The  exact  locations  and  depths  of  the  web 
attachment  fastener  holes  vary  and  in  some  wings  B16  ne^ly  intersects  with  C7. 
For  the  purpose  of  this  analysis  the  centre  lines  of  holes  C7,  A7,  B16  and  Dl5  were 
modelled  as  being  coplanar,  with  B16  intersecting  C7.  The  model  dimensions  are 
shown  in  Figure  2. 

As  well  as  the  uncracked  structure,  two  crack  cases  were  analysed.  In  both  cases  a 
pair  of  cracks  which  initiated  at  the  intersection  of  holes  B16  and  C7  and  propagated 
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above  and  below  B16  were  modelled.  For  case  i  ihe  crack  length  is  lOinin  with 
circular  crack  fronts  above  and  below  hole  U16.  For  case  2  the  crack  has  propagated 
along  the  full  length  of  B16  (13.5mm)  with  a  circular  crack  front  below  and  a 
horizontal  crack  tip  at  the  base  of  the  spar  cap  flange.  The  crack  geometries  are 
shown  in  Figure  3. 

The  analysis  was  performed  using  the  PAFEC  level  6.2  finite  element  analysis  pro¬ 
gram  and  run  on  a  Hewlett-Packard  9000  model  750  workstation.  Three  dimensional, 
twenty-noded  isoparamietric  brick  elements  and  fifteen-noded  isoparametric  wedge 
elements  were  used.  The  aluminium  jdloy  spar  cap  and  steel  wing  attachment  fitting 
were  modelled  separately.  The  bolted  joint  was  approximated  by  modelling  the  steel 
WAF  bolts  so  that  the  bolt  mesh  shared  nodes  with  the  inboard  side  c  ♦he  spar  cap 
holes  and  the  outboard  side  of  the  WAF  holes.  The  F.E.  mesh  is  shown  m  Figure  4. 
It  containes  2764  elements  and  37980  degrees  of  freedom.  A  linear  elcistic  analysis 
was  performed  in  each  case. 

Nodes  in  the  x-y  plane  were  restrained  in  the  z-direction  at  Y884.  This  included  both 
the  spar  cap  and  WAF  and  is  not  strictly  consistent  with  the  actual  deformation  in  a 
bolted  joint.  The  effect  of  this  restraint  was  to  underestimate  the  load  transfer  into 
the  WAF  resulting  in  higher  spar  cap  stresses.  This  is  conservative  for  this  analysis. 
Despite  this,  the  modelling  of  the  WAF  and  bolts  should  produce  a  realistic  stress 
concentration  in  the  WAF  bolt  holes.  Nodes  in  the  z-x  plane  at  y=0  were  restrained 
in  the  y-direction,  see  Figure  4.  This  provides  a  bending  restraint  which  is  consistent 
with  a  full  scale  fatigue  test  [3]  in  which  equal  strains  were  me2isured  at  the  top  and 
bottom  of  the  spar  cap. 

A  uniform  tensile  pressure  load  of  14.5  MPa  was  applied  to  the  spar  cap  cross-section 
at  Y920,  see  [3,4].  This  corresponds  to  the  -fig  normal  load  factor. 

The  spar  cap  material  is  A1  7075,  with  E  =  71000MFa,  v  =  0.33,  Ftu  =  530MPa, 
Fty  =  AbAMPa  and  Kjc  =  Pay/m.  The  WAF  and  bolts  are  steel  with 

E  =  210000MPa  and  v  =  0.3. 


3.  CALCULATION  OF  STRESS  INTENSITY  FACTORS 

The  stress  intensity  factor  was  calculated  at  various  locations  around  the  crack  front 
from  the  FE  displacement  solutions  using  equation  1. 


Ki  = 


uE 

4(1  —  1/2) 


(1) 


where  u  is  the  mode  1  nodal  displacement  near  the  crack  tip,  E  is  the  elastic  modulus, 
V  is  Poisson’s  ratio  and  /  is  the  distance  from  the  node  to  the  crack  tip  [5]. 
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4. 


RESULTS 


4.1  Uncracked  Model 

A  section  of  the  F.E.  mesh  through  the  intersecting  holes  for  the  uncracked  case  is 
shown  in  Figure  5.  For  this  section  the  net  section  stress  is  equal  to  20.2  MPa/g.  The 
maximum  principal  stress  solution  is  shown  in  Figure  6.  The  maximum  principcil 
stress  in  hole  D15  is  52  MPa  which  corresponds  to  a  net  section  stress  concentration 
factor  of  2.6.  For  WAF  bolt  hole  A7  the  maximum  stress  is  39  MPa  giving  a  stress 
concentration  factor  of  1.9,  and  for  the  intersection  of  B16  and  C7  the  maximum 
stress  is  73  MPa  which  corresponds  to  a  stress  concentration  factor  of  3.5.  These 
results  are  summarised  in  Table  1. 

The  residual  strength  of  the  ligament  above  D15  may  be  estimated  by  calculating 
the  point  at  which  the  von  Mises  equivalent  stress,  (<Te),  reaches  the  ultimate  tensile 
stress  of  the  material.  However,  a  significant  amount  of  ductile  yielding  would  be 
expected  at  the  stress  concentration  prior  to  actual  failure.  A  linear  elastic  analysis 
is  therefore  conservative.  For  the  uncraw^ked  case  the  peak  von  Mises  stress  in  the 
ligament  is  40  MPa/g.  The  average  von  Mises  stress  in  the  ligament  is  22.4  MPa. 
The  residual  strength,  Ftulif^elg),  of  the  ligament  is  13.2g. 

The  maximum  von  Mises  stress  at  the  intersection  of  B16  and  C7  is  60MPa/g  which 
corresponds  to  a  residual  strength  of  8.8g. 


4.2  Crack  Case  1 

The  solution  for  maximum  principal  stress  for  crack  case  1  is  shown  in  Figure  7. 
Note  that  the  stresses  for  the  elements  adjacent  to  the  crack  tip  are  not  shown.  The 
di8pljw:ed  mesh  (highly  exaggerated)  is  shown  in  Figure  8. 

The  maximum  von  Mises  stress  in  the  ligament  above  hole  D15  has  increased  to 
60.8  MPa/g,  giving  a  ligament  residual  strength,  FtuUffelg),  of  8.7g.  The  average 
ligament  von  Mises  stress  is  26.6  MPa/g. 

For  the  crack  front  above  hole  B16,  the  stress  intensity  factor  varies  from  a  maximum 
at  the  top  of  3.69  M Pa-^/m/g  to  a  minimum  of  2.48  M Pay/mlg  near  the  centre, 
see  Figure  9.  The  average  value  of  Kj  for  the  upper  crack  is  2.8  MPay/in/g.  For 
the  crack  front  below  B16  the  stress  intensity  factor  varies  from  a  maximum  of 
3.37  MPoy/mfg  at  the  top  to  2.1  MPay/m/g  at  the  bottom.  The  average  is  2.3 
MPoy/mfg  .  The  residu2d  strength,  Kjcl{Kj  jg),  based  on  average  and  maximum 
stress  intensity  factors  are  11. 6g  and  8.8g  respectively. 
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4.3 


Crack  Case  2 


A  section  through  the  F.E.  mesh  for  crack  case  2  is  shown  in  Figure  10.  The  solution 
for  maximum  principal  stress  is  shown  in  Figure  1 1  and  the  displaced  mesh  is  shown 
in  Figure  12.  The  maximum  von  Mises  stress  in  the  ligament  above  hole  D15  is  78.7 
MPa/g,  giving  a  residual  strength,  Ftu/{<^el9),  of  6.7g.  The  average  ligament  von 
Mises  stress  is  32.3  MPa/g. 

For  the  upper  crack  front  the  stress  intensity  factor  varies  from  a  maximum  of  3.38 
A/ Pay/mfg  in  the  centre  to  3.10  M Pay/mfg  at  the  ends.  The  average  value  is  3.25 
M Poy/m/g.  In  the  lower  crack  front  the  stress  intensity  varies  from  a  maximum  of 
4.9  M Pay/mlg  at  the  top  edge  to  a  minimum  of  2.52  M Pay/mfg  at  the  centre.  The 
average  stress  intensity  is  3.06  M Pay/mfg,  see  Figure  13.  The  residual  strength, 
f g),  for  average  and  maximum  stress  intensity  factors  are  10. Og  and  6.7g 
respectively. 

The  results  for  crack  cases  1  and  2  are  presented  in  Table  2. 


5  PLASTIC  ANALYSIS 

The  linear  elastic  analysis  predicted  a  peak  stress  at  hole  D15  much  greater  than 
the  average  ligament  stress,  indicating  that  the  stress  concentration  is  very  localised. 
A  significant  amount  of  plastic  yielding  and  associated  stress  redistribution  would 
occur  before  tensile  failure  of  the  ligament.  To  predict  accurately  the  point  at  which 
the  von  Mises  stress  reaches  the  material  ultimate  tensile  stress  it  is  necessary  to 
account  for  this  material  non-line^lrity. 

A  plastic  aneilysis  of  the  structure  around  hole  D15  was  therefore  performed.  For 
each  of  the  linear  elastic  models  of  the  spar  cap,  the  elements  in  the  spar  cap 
flange  around  hole  D15  were  extracted  and  the  mesh  around  the  hole  refined,  see 
Figure  14.  The  nodal  displacements  from  the  linear  elastic  solutions  were  applied 
around  the  boundary.  This  modelling  strategy  is  valid  provided  the  influence  of  the 
localised  yielding  is  small  at  the  model  boundary.  The  tensile  stress-strain  curve  for 
aluminium  alloy  7075  [6]  was  approximated  by  a  series  of  straight  lines  as  shown  in 
Figure  15. 

The  non-linear  solution  was  run  using  0.5g  increments  and  the  load  required  to  reach 
material  failure  determined.  This  was  defined  ais  the  point  at  which  the  maximum 
von  Mises  stress  equalled  the  material  ultimate  tensile  stress. 

The  von  Mises  stress  solution  for  the  uncracked  case  is  shown  in  Figure  16.  The 
plastic  region  is  small,  extending  over  about  20  percent  of  the  ligament.  The  use  of 
linear  elastic  displacements  at  the  boundary  is  therefore  a  reasonable  approximation. 

The  ligament  residual  strength  results  are  listed  in  Table  3.  Three  definitions  of  fail¬ 
ure  are  compared.  In  the  first  definition,  failure  occurs  when  the  average  ligament 
stress  equals  the  material  yield  stress.  In  the  other  two  definitions  the  ligament 


is  deemed  to  have  failed  when  the  maximum  von  Mises  stress  reaches  the  mate¬ 
rial  ultimate  stress,  one  for  the  linear  elastic  solution  and  the  other  for  the  plastic 
solution. 


6.  DISCUSSION 

For  the  uncracked  structure  the  highest  stress  lies  at  the  intersection  of  holes  B16 
and  C7  and  corresponds  to  a  stress  concentration  factor  of  3.5.  The  intersection 
point  is  therefore  a  possible  location  for  crack  initiation. 

For  crack  case  I  the  upper  crack  front  geometry  appears  to  be  realistic  since  the 
values  of  stress  intensity  are  similar  at  each  end.  The  variation  of  A'/  around  the 
crack  front  is  consistent  with  the  solution  for  a  quarter-circular  crack  in  a  square 
bar,  see  [7].  For  the  lower  crack  front  the  value  of  Kj  at  the  top  is  significantly 
higher  than  at  the  bottom.  This  is  due  to  the  stress  concentration  at  hole  B16. 
Crack  growth  is  therefore  likely  to  propagate  f2ister  along  B16  and  result  in  a  more 
elliptical  crack  front.  The  maximum  value  of  Kj  reaches  the  critical  value  at  8.8g 
(lower  crack)  and  the  average  A'/  becomes  critical  at  11. fig  (upper  crack). 

For  crack  case  2  the  maximum  stress  intensity  factors  have  increased  to  give  a 
residual  strength  of  lO.Og  (lower  crack)  and  the  average  value  results  in  a  residual 
strength  of  fi.7g  (upper  crack). 

Residual  strength  verses  crack  length  curves  for  fracture  of  the  proposed  cracks  and 
for  tensile  failure  of  the  ligament  above  D15  are  shown  in  Figure  17.  There  was  a 
significant  variation  of  the  stress  intensity  factor  along  the  crack  fronts,  so  results 
based  on  both  the  maximum  and  average  stress  intensity  factors  are  included.  Only 
the  ligament  residual  strength  based  on  the  plastic  analysis  is  included  since  the 
linear  elastic  result  was  considered  too  conservative  for  the  purpose  of  this  analysis. 

The  residual  strength  of  the  spar  cap,  based  on  the  maximum  stress  intensity  factor, 
is  below  the  residual  strength  of  the  ligament  and  falls  below  the  limit  load  of  8g  for  a 
crack  length  of  1 1.3  mm.  For  comparison,  the  residual  strength  based  on  the  average 
stress  intensity  factor  lies  above  the  ligament  residual  strength.  There  is,  therefore, 
some  uncertainty  in  predicting  the  load  at  which  fr2M:ture  would  occur.  However, 
these  results  do  indicate  that  failure  could  occur  from  the  cracked  region  in  the 
presence  of  an  intact  flange  hole  ligament:  in-service  NDI  monitoring  of  the  flange 
surface  above  the  ligament  cannot,  therefore,  guarantee  the  structural  integrity  of 
the  section. 


7.  CONCLUSIONS 

1.  The  analysis  indicates  that  failure  of  a  cracked  section  could  occur  with  an 
intact  ligament. 

2.  NDI  of  the  spar  cap  flange  cannot  guarantee  the  structural  integrity  of  the 
section. 
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Table  1. 


Stress  concentration  factors,  K,.  for  the  uncracked  structure 


Location 

Max.  Principal  Stress 
(MPa) 

Kt 

Upper  web  fastener  hole 
D15 

52 

2.6 

WAF  bolt  hole 

A7 

39 

1.9 

Intersection  of  holes 
B16andC7 

70 

3.5 

Table  2.  Stress  Intensity  Factors,  K|  (MPaVm  /g) 


Crack  case 

Upper  crack 

Lower  crack 

Maximum 

Average 

Maximum 

Average 

1 

3.7 

2.8 

3.4 

2.3 

2 

3.4 

3.3 

4.9 

3.1 

Table  3.  Ligament  residual  strength  (g) 


Average  Ligament 

Material  Ultimate 

Material  Ultimate 

Case 

Yield 

Stress  -  Linear 

Stress  -  Plastic 

uncracked 

20.1 

13.3 

17.0 

crack  case  1 

16.9 

8.7 

10.6 

crack  case  2 

13.9 

6.7 

8.8 

i 


LOWER  SPAR  CAP  DEUTSCH  FASTENERS  ^ 

WOT  TO  SC*Lt 


Figure  1  Lower  spar  cap  and  WAF  assembly 


Figure  2 


F.E.  Model  Dimensions 


Figure  8  Displaced  mesh  for  crack  case 


Figure  9  Stress  intensity  factors  for  case  1 


Figure  10  F.E,  mesh  for  crack  case  2 


Figure  11  Maximum  principal  stress  solution  for  crack  case 


Figure  12  Displaced  mesh  for  crack  case  2 
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K|  (MPaVm/g) 
3.10 
3.38 
3.14 
4.90 
2.85 
2.52 
3A1 


Location 

1 

2 

3 

4 

5 

6 
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Figure  13  Stress  intensity  factors  for  case  2. 


Plastic  solution  (von  Mises  stress)  of  the  flange  hole  for  the  uncracked  case 
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